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Evolution of Fragments Formed at the Rupture of a Knotted Alkane Molecule 
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Abstract: Common experience tells us that a knot significantly weakens the polymer strand in 
which it is tied, which in turn leads to more facile chain rupture under tensile loading. Using first- 
principles molecular dynamics calculations we describe the dynamical evolution of the radicals that 
form after chain rupture of a single knotted alkane molecule in their very early stages of life. They 
are able to recombine, to form cyclic alkanes and to undergo disproportionation phenomena with 
nearby chain segments. The breaking of a single knotted polymer chain under mechanical loading 
is thus predicted to reveal phenomena falling in the domain of ultrafast spectroscopy. 



C^ 



X3 

o 



> 

O 

On 
On 
-I— > 



I 

O 

o 



X 



1. Introduction 

The topological properties of knots ||l|, their occur- 
rence in organic and bio-polymers H, and the modifi- 
cations to physical and chemical properties induced by 
the presence of knots Q, are well-established fields of 
research M. Simple single-chain knots are consistently 
observed in real polymers 0, including DNA H and 
other bio-polymers Jq] , and are widely studied with com- 
puter simulations M. In particular, the "trefoil" knot, 
familiarly referred to as the "overhand knot" , is the sim- 
plest and most commonly occurring kind of single-rope 
knot S. Even so, relatively little is known about knots 
in polymer strands at the atomic level M. The response 
of knotted polymers to mechanical deformations, such 
as stretching (tensile loading), has a major impact on 
polymer technology. In fact, techniques, such as atomic- 
force microscopy (AFM) ||l^ or optical tweezers [|ll| , can 
now manipulate single molecules and thus apply tensile 
load at a molecular level. Recent first-principles quantum 
mechanical calculations ||l2[ have shown that the prop- 
erties of a "trefoil" knot in a polyethylene chain, sub- 
jected to tensile loading, closely follow the well-known 
behavior of macroscopic ropes |g]. Indeed, analysis of 
the strain energy distribution along the knotted poly- 
mer chain in terms of a classical force field model p3[ , 
shows that the bonds at the entrance and exit points are 
far more stressed than the others, and that the bonds 
near the center of the knot are still reasonably close to 
their ground state geometry. First principles calculations 
also confirm our common experience that the presence of 
the knot significantly reduces the tensile strain that can 
be stored in the strand before breaking. Moreover, the 
break point in a microscopic polymer strand occurs at 
either the entrance or exit point of the knot |l2| , exactly 
as in a macroscopic rope B , and as recently observed in 
optical tweezers experiments [0 . 

The focus in the present article, is not directed at the 
behavior of the stress distribution and the nature of the 
break point in a knotted polymer but rather on what 



happens to the resulting chain fragments after the break 
occurs. Since the rupture of a polymer strand necessarily 
generates fragments that are radicals one would antici- 
pate interesting chemical phenomena e.g., radical recom- 
bination, which have no obvious corresponding counter- 
part with a macroscopic rope. Accordingly, we have in- 
vestigated the breaking of an individual alkane molecule 
that contains a simple trefoil knot focusing on the dynam- 
ical evolution of the radicals that form after the break. 
Although this is an idealized system, studying the nature 
of chemical processes characterizing individual molecules 
is a necessary first step in understanding polymer prop- 
erties at the macroscopic level. 

2. Computational details 




Bond number 

FIG. 1. Strain energy distribution for a C28 n-alkane 
molecule that has been tied into a trefoil knot. The actual 
configuration, which has an end-to-end separation L — 14. OA, 
is shown in Fig. 2a. This representative configuration was ob- 
tained from a classical MD simulation using a force field model 
(8). The system was evolved for 200 ps and then annealed for 
50 ps. 

We have shown elsewhere ||l2| that an n-alkane con- 
taining 28 carbon atoms and, accordingly, 58 hydrogens is 



large enough to be knotted into a trefoil. In that work, as 
in the present case, tension was applied by constraining 
the end-to-end length of the molecule, i.e. the Ci — C28 
distance L to a series of fixed values. Calculations with 
gradually increasing L were performed until chain rup- 
ture occurred. Classical molecular dynamics (MD) simu- 
lations were used to pre-equilibrate the system for about 
250 ps at any new value of L, thus bridging any two suc- 
cessive quantum-mechanical calculations. The minimum 
end-to-end distance between the terminal methyl group 
carbon atoms L required to produce chain rupture on 
the timescale of the calculations in the knotted C28H5S 
alkane was 13.5 A. This corresponds to a contour length 
of the chain of ca. 41 A, as compared to an equilibrium 
chain length of ca. 35 A. 

The knotted polymer chain was studied using Car- 
Parrinello [Q molecular dynamics |1^ , based on density- 
functional theory (DFT) with Becke and Lee- Yang-Parr 
(BLYP) [Q for the exchange and correlation functionals. 
A valence electron pseudopotential/plane-wave scheme 
was employed, with a kinetic-energy cutoff of 60 Ry, a 
F-point sampling of the Brillouin zone. Since chain rup- 
ture yields unpaired electrons, the spin variables were 
handled through the so-called local spin-density approx- 
imation |l^ . Periodic boundary conditions were used to 
the aim of extending the present single-chain results to 
liquid-crystal systems at the same level of accuracy and 
without any loss of consistency. 

First-principles MD is particularly useful in the study 
of molecular processes involving breaking and forma- 
tion of chemical bonds. Even though DFT-based esti- 
mates of transition states energetics have only limited 
accuracy [|l8[ , the qualitative estimates of relative values 
of barriers between different reaction paths is generally 
more satisfactory 19 1. In addition, we have shown else- 
where Q that the methodology employed yields bond 
energies accurate to better than 2 %. In the present case, 
as it will be shown in the following, the chemistry of gas- 
phase-like chain radicals in their earliest stages of life is 
dominated by barrierless reactions such as recombination 
and disproportionation. Even though the latter process 
is energetically less favorable than the former, the ratio of 
the respective reaction rate constants is about 0.4-0.6 for 
ethyl radicals, and even larger for longer fragments. This 
is due to an activation barrier that, for both reactions, 
is zero within experimental error [E0[ . The occurrence of 
one of them, rather than the other, appears to be purely 
driven by dynamical factors, such as steric effects and 
topological constraints. The CPMD approach allows one 
to follow dynamical evolution of systems with up to a few 
hundred atoms. Interatomic forces are determined by the 
instantaneous electronic structure. Unlike classical force 
fields, this feature ensures an accurate description of very 
different situations such as, in the present case, changes 
in hybridization status and coordination |gl|. To fol- 
low the dynamical evolution of the system at room tem- 



perature (T — 300K), we employed a time step of 6.0 
a.u. {— 0.145 fs). The starting configuration employed 
in the CPMD calculations, was obtained from a series 
of classical MD simulations with an empirical force field 
model ijl^. We discuss in detail the situation shown in 
Fig. n^, which corresponds to a large tensile load at the 
entrance and exit points of the trefoil 1 22 1 . 



3. Results 

3.1 Double chain rupture and cycloalkane for- 
mation. The corresponding configuration is shown in 
panel A of Fig. g. Panels B-F depict the initial dynami- 
cal evolution of the carbon backbone of the knotted poly- 
mer. In less than 100 fs, the long chain molecule breaks 
at both the entrance and exit of the trefoil (panels B and 
C), thereby generating a large diradical b-c (following 
the labeling scheme of panel C) and two smaller radicals 
a and d. The system is simply unable to redistribute 
the imposed tensile stress (see Fig. ^) quickly enough 
to avoid a second bond dissociation Q. Slightly more 
than half of the total strain energy (w 300kcal/inol) is 
released into radical formation energy (2 * 8Skcal/mol). 
The remaining part is transformed into kinetic energy of 
the radicals, that strongly recoil after the chain actually 
"snaps" . 







t=160fs 



t=200 fs 




t=260 fs Y X 



FIG. 2. The snapshots show the first 260 fs of dynam- 
ical evolution of the system. Only the carbon atom back- 
bone is drawn. Panel A gives the definition of L as the fixed 
end-to-end distance. In panels B-E the "interesting" carbon 
atoms are enlightened, and labeled with letters in panel C 
(see text). The axis orientation is given in panel F. 

After the generation of the three radicals, one might 
have expected to observe the radical recombinations a- 



c and b-d, thus creating an unknotted strand from 
the knotted one. However, instead, the radicals are so 
strongly reactive that, within about 150 fs after the 
chain rupture, the carbon atoms labeled as b and c form a 
chemical bond (panel E) . This phenomenon is favored by 
the previously mentioned recoiling of the b and c frag- 
ments after the rupture, that come very close to each 
other and thus recombine. This reaction results in the 
formation of a CigH^s cyclic alkane, still leaving two 
active radical fragments, one C^Hj- (a) and the other 
CqHi^- (d), whose terminal methyls remain anchored to 
their initial positions, separated by L = 14.0 A. The 
"constraint" on the terminal methyls, thus inhibits the 
direct a-d recombination (panel F). This is precisely the 
situation that obtains if a knotted polymer is broken us- 
ing a laser tweezers technique ^J. 

3.2 Disproportionation. In the case of an isolated 
knot, in order to have a reaction between the fragments 
a and d, it will be necessary to remove the constraint on 
the terminal methyls. However, due to periodic boundary 
conditions, the simulation system we employ is not com- 
pletely isolated, it is surrounded by its periodic images. 
In the present case, the box length in the z direction, i.e. 
the direction along which the linear ends of the knotted 
chain are aligned, is short enough to allow the two bro- 
ken radicals to interact with their periodic images after 
the knot breaks. This happens because, after the trefoil 
ruptures, the fragments recoil strongly along the z direc- 
tion and the —z direction, respectively. In particular, as 
shown in panels C-F, the z coordinate of the atom la- 
beled as a tends to increase, while the opposite is true 
for the d atom. 

Fig. y illustrates the subsequent dynamical evolution 
of the system, starting about 700 fs after the ring closure 
that generated the cyclic alkane. It is noteworthy that 
the longest of the two fragments (whose unsaturated car- 
bon atom is labeled as d) is actually the periodic image, 
in the z direction, of the bottom radical of Fig. g (pan- 
els C-F). A simple reaction that was expected to occur 
was the saturation of the dangling bonds through the re- 
combination of the a carbon with the periodic image of 
the d atom. Instead, the shorter radical attacks the car- 
bon bonded to the d atom (panel B), rather than the 
d carbon itself, "stealing" one of its hydrogens (panels 
C-D) and then recoiling away (panels E-F), in a typical 
disproportionation reaction. 

This reaction takes place in a very short time 
(w 50 fs), leaving the two terminal atoms of the longer 
fragment (displayed in dark gray) to form a double bond, 
thus saturating all the unpaired electrons. In panel F the 
sp^ hybridization status of these two carbons is clearly 
observable from their coplanarity with the hydrogens to 
which they are bonded. The energy gain due to the dis- 
proportionation is ca. 55 kcal/mol, which compares fa- 
vorably with the experimental enthalpy of the analogous 
reaction between two ethyl radicals (60 kcal/mol). Due 



to the high computational demand of first-principles cal- 
culations, we could only follow a few other trajectories. 
However, in each case we observed the same phenomena. 
Even though the results may not be meaningful in a strict 
statistical sense, the energetics of the activation barriers 
of the competing processes suggest that, on an ultrafast 
timescale, the occurrence of such reactions is dominated 
by steric effects, which, in the present case, tend to favor 
disproportionation against recombination. As expected, 
the cyclic alkane ring is stable and does not take part in 
any of these subsequent reactions. 
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FIG. 3. Snapshots of the dynamical evolution of the sys- 
tem after the formation of the cyclic alkane. Carbon atoms 
are displayed in gray and hydrogens in white. One of the two 
active carbons and the hydrogen atom that is to be abstracted 
are shown in black. The other unsaturated carbon, and its 
bonded C atom are displayed in dark gray. In panel A the 
two active carbons are labeled according to the notation of 
Fig. 2c, and the axis orientation is given. 

The nature of the relevant bonds can be actually mon- 
itored throughout the simulation. In Fig. we show two 
meaningful atom-atom distances as function of time. In 
panel A the bond length of the two carbon atoms col- 
ored in dark gray in Fig. 2 is displayed. The three 
distinct situations that the system experiences can be 
clearly identified in this picture. In the very first part of 
the dynamical evolution (« 100 fs), the C — C distance 
is larger than the single bond equilibrium length (1.54 
A) . Here, although all the carbon atoms of the chain are 
still bonded, the whole system is highly distorted and all 
the bonds quite stretched. After ca. 200 fs the knot 
breaks, and the C — C distance then oscillates around an 



average value of 1.49 A, which is the typical bond length 
between a sp^- and a sp^ -hybridized carbons. Shortly 
after 1100 fs of dynamics, the "black" hydrogen of Fig. 
2 is abstracted by the unsaturated carbon atom of the 
shorter radical (panel B), forming a stable bond with a 
large-amplitude vibration that relaxes in a time of about 
200 fs. In this case, part of the rototranslational kinetic 
energy of the radicals is converted into molecular vibra- 
tions of the shorter a fragment, mainly the C-H bonds 
of the methyl group involved in the reaction. 

Concomitant with the disproportionation reaction, the 
"dark gray" C — C bond length shortens and begins to 
oscillate around an average distance of 1.34 A and at a 
higher frequency. This behavior, observable in the right 
part of panel A of the picture, is characterized by a 
shorter-period vibration. 
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FIG. 4. Panel A: time evolution of the bond length of the 
atoms highlighted in dark gray in Fig. 3. Panel B: time evo- 
lution of the distance between the carbon and the hydrogen 
atoms displayed in black in Fig. 3. The inset (panel C) shows 
an enlargement of the large-amplitude C — H vibrations fol- 
lowing the abstraction reaction. 

4. Summary 

In conclusion, we report herein a first-principles MD 
study of the dynamical evolution of hydrocarbon radicals 
generated by the breaking of a single polyethylene-like 
chain containing a simple trefoil knot. The present work 
suggests that the early stages of evolution are likely to 



generate interesting chemistry, due to the resulting un- 
saturated alkane fragments. The 1.5 ps time evolution 
we monitor is characterized by barrierless reactions that 
are biased by steric effects. The velocity distribution of 
the recoiling radicals at the rupture is a key factor, and 
is intrinsically related to the topological constraints im- 
posed by the presence of a trefoil knot. The most logi- 
cal chemical reaction, namely chain recombination, is to- 
tally bypassed in favor of ultrafast phenomena including 
the formation of a diradical |2^ that generates a cyclic 
alkane, followed by disproportionation, and simultane- 
ous carbon-carbon double bond formation. Observation 
of the phenomena reported herein is likely to present a 
challenge for experimentalists, particularly because the 
reactions occur in the femtosecond range [g6|. The re- 
sults we currently report do not include the effects of 
surrounding neighboring chains on the chemistry of the 
fragments. However, preliminary calculations on a knot- 
ted alkane embedded in a nematic liquid crystal of like 
molecules show that the environment does not seem to 
significantly modify the observed chemical reactions. 
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